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Abstract 



Using a Dirac-matrix substitution rule, applied to the electric charge, the anomalous mag- 
netic moments of fermions are incorporated in local form in the two-body relativistic wave 
equations of constraint theory. The structure of the resulting potential is entirely deter- 
mined, up to magnetic type form factors, from that of the initial potential descibing the 
mutual interaction in the absence of anomalous magnetic moments. The wave equations 
are reduced to a single eigenvalue equation in the sectors of pseudoscalar and scalar states 
(J = 0). The requirement of a smooth introduction of the anomalous magnetic moments 
imposes restrictions on the behavior of the form factors near the origin, in x-space. The 
nonrelativistic limit of the eigenvalue equation is also studied. 
PACS numbers: 03.65.Pm, 11. 10. St, 12.20.Ds. 
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1 Introduction 



The two-body relativistic wave equations of constraint theory have the main feature of 
describing the internal dynamics of the system by means of a manifestly covariant three- 
dimensional formalism |], 0, |3|, f|, ||, ||, 0, @, 0; relative energy and relative time vari- 
ables are eliminated there through constraint equations. These wave equations, which 
can be constructed from general principles, have also the property of allowing a three- 
dimensional reduction of the Bethe-Salpeter equation by means of a Lippmann-Schwinger- 
quasipotential type equation fllUL O, O, O, [T3, [15], [TBI, IT7L |18]] that relates the two-body 



potential to the scattering amplitude [^, |20[ . In this way, the potential becomes calcula- 
ble, in perturbation theory, from Feynman diagrams. 

In a recent work |?TJ , we applied this calculational method to the evaluation, in certain 
approximations, of the potentials in the cases of scalar and vector interactions, mediated 
by massless photons. It turns out that at each formal order of perturbation theory, 
which is now reorganized by the presence of additional three-dimensional diagrams due 
to the constraints, the leading infra-red terms are free of spurious singularities and can be 
represented in three-dimensional x-space as local functions of r, proportional to (g 2 /r) n , 
where r is the cm. relative distance, g the coupling constant and n the formal order 
of perturbation theory. The series of leading terms can be summed and result in local 
functions (in r) for the expressions of the potentials. The latter are compatible with the 
potentials proposed by Todorov in the quasipotential approach on the basis of minimal 
substitution rules |1| and later investigated in the fermionic case by Crater and Van 
Alstine 0. 

The purpose of the present paper is to take into account, in the case of vector interac- 
tions, the effects of the anomalous magnetic moments of fermions. In usual calculations 
in QED, because of the smallness of the coupling constant, the latter are evaluated at 
leading order of perturbation theory in the nonrelativistic limit (0(a 5 ) effect, where a 
is the fine structure constant). However, in strong coupling problems, like in the strong 
coupling regime of QED or in the short- distance (vector) interactions of QCD, nonper- 
turbative contributions of the anomalous magnetic moments may become sizable and the 
incorporation of higher order effects becomes necessary. 

In order to include the main effects of the anomalous magnetic moments into the 
potentials, we evaluate their contributions through the vertex corrections. In the lowest 
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order graph the anomalous magnetic moment appears by means of a substitution rule that 
replaces each charge (coupling constant) by a Dirac-matrix function. We then introduce 
this typical vertex correction at the vertices in each order of the perturbation series of 
the vector potential determined previously [plfl . Although the substitution rule utilized 
above is rather simple, it generates in the higher order terms technical complications for 
the summation of the perturbation series of the potential. The reason for this is the 
new Dirac-matrix structure that results from the higher order terms. Up to now, all the 
potentials that were considered in the constraint theory wave equations had dependences 
on the Dirac matrices only through pairs of 71 and 72 matrices, the indices 1 and 2 referring 
to the two fermions, respectively (like 71.72, 715725, etc.)— a feature that considerably 
simplifies many algebraic operations as well as the reduction process to the final eigenvalue 
equation. The presence of the vertex corrections, even though globally symmetric in the 
exchanges 1 <-> 2, breaks this symmetry in the individual terms and introduces new types 
of structure not present in previous calculations. It is the presence of these terms that 
makes calculations rather complicated. While the potential can still be represented in 
a somehow compact form, the final eigenvalue equation for general quantum numbers 
becomes less easy to handle. It takes a relatively simple form only in the sectors of j = 
states (pseudoscalar and scalar), to which we have limited our final analysis. The ground 
states of these sectors are precisely those which may be concerned with spontaneous 
breakdowns of symmetries (chiral and dilatational). 

The plan of the paper is as follows. In Sec. 2, we consider lowest order perturbation 
theory and determine the substitution rule to be used for the vertex correction. In Sec. 
3, the vertex corrections are incorporated into the vector potential. The new form of 
the latter is determined in Sec. 4, by resumming the corresponding perturbation series. 
In Sec. 5, the wave equations are reduced, for the j = states, to a final eigenvalue 
equation. In Sec. 6, we analyze the effects of the anomalous magnetic moments in several 
limiting cases. The requirement from the accompanying form factors of not aggravating 
the singularities of the initial potential, leads to restrictions on their behavior near the 
origin in x-space. The nonrelativistic as well as the one-particle limits are also checked. 
Concluding remarks follow in Sec. 7. 
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2 Structure of the two-body potential with 
anomalous magnetic moments 



In order to determine the structure of the potential in the presence of anomalous magnetic 
moments, we start from the expression of the coupling of a pointlike particle (fermion 
1) with charge e\ and anomalous magnetic moment K\ to an external electromagnetic 
potential and to its field strength tensor F^ v : 

(Ti^ + ^i^"") (2-1) 

\ a iiv — 2?[7/f 7"]]- In the case of a mutual interaction with another particle 2, expression 
(|2.1| ) represents the lowest order perturbation theory result, where potential is itself 
expressed in terms of the photon propagator and its coupling to particle 2 (which we 
suppose for the moment without anomalous magnetic moment): 

A, = D^Y 2 , (2.2) 

where D^ v is the photon propagator including the coupling constants at its two ends. 
In the three-dimensional formalism of constraint theory, the Bethe-Salpeter kernel is 



projected on the constraint hypersurface and the wave function expanded around it [19], 
pOfl . in the cm. frame, this amounts to projecting the kernel on the hypersurface where 
the temporal component of the momentum transfer is zero. In a covariant formalism, one 
first decomposes four-vectors along transverse and longitudinal components with respect 
to the total momentum P: 

P = Pi + P2 , P = g(pi - P2) , X = ^(x 1 + x 2 ) , X = Xi - x 2 , 

xl = X M - P.XP^ , X L = P.X , P L = , P M = , 



x = x — (P.x) , r = y — x T2 , 

TJ = 7m - P-lh , 1l = P-l , M = mi + m 2 . (2.3) 

Thus, in constraint theory, the propagator in Eq. (|2.2|) depends on x, in x-space, through 
x T only; in the Feynman gauge, to which we stick throughout this work, it has the 
expression (in lowest order) [ F2"T| : 



= g^D(x T2 ,P L ) , D{x T \P L ) = ■ (2.4) 



Using Eqs. ( |2.2| ) and Q2.4| ) the field strength tensor takes the form: 



d^A u — dyA^ 



where the dot operation represents derivation with respect to x T2 \ 





(2.5) 



/ = 



Q X T2 



f 



[2.6) 



Expression (|2.1[ ) then becomes: 



1 K\ 

7im^ + 2 Kl(Tl ^ F>1V = D t7i-72 - ^(7i T -^ T 7i-72 - 7i-727?^ T )] 

= ~D(1 - i— 7l T . — )7 1 .7a 
2 r r 

1 n/i ' Kl t ^ l f 
-D(l - . — J71.72 710720 

2 r r J 



+7io720 



(2.7) 



[The dagger represents hermitian conjugation.] We deduce that to this order the anoma- 
lous magnetic moment appears through the following matrix substitution of the charge 
ex. 

e, - e[ = ei (l-t^l^) . (2.8) 

The above calculations can be repeated at the particle 2 (antifermion) vertex, where the 
charge substitution becomes: 



( D 2 



K2 T X 



e' 2 = e 2 (l -i— 7 2 



(2.9) 



[ea is the fermion 2 charge; the passage to the antifermion is obtained in momentum space 
with the replacement p\ — > —p 2 , which yields for the momentum transfer q = pi — p[ — > 
— p 2 +P2 = 9 an d hence in x-space x — > x; the Dirac matrices 72 act on the wave function 
from the right.] 

The mutual interaction potential then becomes to lowest order: 



V = \d{1 - - ^ 72 T -)7i-72 

2 r r r r 



+710720 



~£>(1 - z— 7f .— )(1 - z— 7 2 '-— )7i-72 
2 r r r r 



. ^2 t 



t 

710720 



(2.10) 



It is natural to generalize the substitution rules ( |2.8|) -( [2T9"D to higher orders, with the 
difference that the lowest order anomalous magnetic term n/r should now be replaced by 
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a form factor b(r): 

1 - z&WtJ-^t) (a = 1,2). (2.11) 



e' 



The form factor 6(r) should ensure for the anomalous magnetic term a smooth behavior at 
the origin, in order not to enhance the existing singularities of the propagator D. Taking 
into account the lowest order expression of k, k = a convenient parametrization for 

b{r) is: 

2mr + (to/(27r))c a (r) 
e 2 

a = — , ei = e 2 = e . (2-12) 

47T 



Expression (|2.10 ) is then generalized to the following form: 



1 / HP HP HP \ 

V = \{A- iB xl \ X - - rB 2l T 2 X - - C 7 f.^ 7 J.^) Tl . 72 

2 V r r r r J 

j / ry> r f~' ry'-f iy*^ ty*~^ \ 

+ i 7l . 72 U + ^ l7 f.^+^ 272 ^-C 7 r^ 72 r .^) , (2-13) 
2 V r r r r J 



where the potentials A, B%, B 2 and C are completely determined, by means of the sub- 
stitution rules ([2.111) ) from the expression of V in the absence of anomalous magnetic 
moments. The latter expression has the form: 

V = 4)71-72 . (2.14) 

[Aq is denoted by V 2 in Ref. |2TJ.] 

For the Todorov potential 0, 0, 13> is: 

4 = + (2.15) 

However, other effective expressions could be used for A Q as well. 

Actually, the potentials that appear in the constraint theory wave equations are func- 
tions of V [Eq. Q2.13| )1 through exponentiations; therefore, we shall need to calculate such 
exponential functions. This is the main content of Sec. 3. 
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3 Wave equations 



The wave equations of constraint theory for a fermion-antifermion system can be written 
in the form 



22 



(ji-Pi - mi) * = (— 72.P2 + m 2 )V \1> 
[-J2-P2 - m 2 ) * = (71 -Pi + rrii)V § , 



(3.1) 



where \I> is a spinor function of rank two, represented as a 4 x 4 matrix function; the Dirac 
matrices (72) of the antifermion act on \1/ from the right; the total and relative variables 
were defined in Eqs. (|2.3j ); V is a Poincare invariant mutual interaction potential. 
Equations ( |3.1|) imply the constraint 



(p? - pI) - ( m i - m l) 







or equivalent ly 



C(p) = 2P lPl 



fm? — mo) ~ 







(3.2) 



(3.3) 



which allows the elimination from the wave equations of the relative longitudinal mo- 
mentum in terms of the masses and the cm. total energy. The wave function for 
eigenfunctions of the total momentum P, has then the structure: 



*{X,x) = e -iP-X e -i{m x -m 2 )x L /{2P L ) 



(3.4) 



The positivity conditions of the norm of \& imply that V should satisfy the inequality 
jTrV^V < 1 |8|, 22fl . A convenient parametrization satisfying this inequality for potentials 
commuting with ^nl2L was proposed by Crater and Van Alstine [p3| ; it is: V = tanh V. 
It turns out that the perturbation series of the leading infra-red terms in QED in the 



Feynman gauge provides a potential V that is compatible with this parametrization [21 



(cf. Eqs. (|2.14j) - (|2.15|) , where V is denoted by Vo). For more general potentials which do 
not commute with 71L72L, the generalization of the above parametrization is [^2]| : 



Iil^lV = tanh(7i L 72L^) • 



Equations ( |3.1| ) are then transformed with the change of wave function: 



(3.5) 



^ = cosh( 7lL 7 2L V) * 



(3.6) 
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they become: 



(71 .p x - mi) cosh(7 1L 7 2L V A ) ^ = (-72^2 + ^)7il72L smh(7 1L 7 2L ^) ^ , 
-12-P2 - m 2 ) cosh(7i L 72L^) * = (71 -Pi + ™i)7ll72L sinh(7i L 7 2 L^) * • (3.7) 



One can also equivalently work in the "Breit representation". Defining 

V B = 1il12lV , 

-v Bx 



ty B = e 



one shows that Eqs. ( |3.1| ) or ( |3~7| ) reduce to thye Breit type equation p2| : 

P L e 2VB - (Hi + ft 2 ) * 



(3.8) 
(3.9) 



(3.10) 



provided constraint (|3T2|)- (|3~3| ) is used; here, Ti\ and TC2 are the covariant free hamiltoni- 
ans: 



H 2 = -m 2 'j2L ~ l2Lll-Pl 



(3.11) 



The normalization conditions of the wave functions \1/ and ^ B were presented in 



Ref. 122 



To solve the wave equations one decomposes the sixteen-component (4 x 4) wave 
function ip along four-component (2 x 2) wave functions: 



1p = 1pl+ ^L^2 + 75^3 + 7^75^4 , 



and similarly for the Breit type wave function ip B : 



IpB = 4>B1 + lL^B2 + 75^B3 + Ihl^Bi ■ 



(3.12) 



(3.13) 



These components are obtained with the projectors [g2 



Vi = ^(1 +7il72l)(1 + 715725) , V 2 = ^(1 +7il72l)(1 -715725) 



V 3 = -(1 -7il72l)(1 + 715725) 



Vi = -(1 - 7il72l)(1 - 715725) • (3.14) 
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The spin operators, which act in the four-component wave function subspaces, are 
defined by means of the Pauli-Lubanski operators: 
h 1 

W 1Sa = —£Z a pp,P P 0%' , W 2Sa = —^aPt/vP^d? ( e 0123 = +1) , 

TI7 TlP-L T TTT ^PL T 

luWlSa = -^-7la7l5 , l2hW 2 S a = ~l^l 2a l2b , 
Wis = = ~\tfP 2 , W S = ^15 + ^ 25 , 

2 \2 4 

w = ' 7~5~ w is-W 2S — > - 72S1.S2 , 
v UPl J c - m - n 

( 2 \ 2 W ls .x T W 2S .x T 4 (si.x)(82.x) 



Wl2 \hP L J x T2 h 2 x 2 

w\ 2 = 1 , 1012(10 - wi 2 ) = w - w 12 ■ (3.15) 



It is clear, from Eqs. (|3.7| ) and ( p . 1 0|) , that one has to calculate the exponential of 



7il72L^, with V having the general structure ( 2.13|) . We have: 
V B = 7il72L^ 

1/ rr>'-l~' nn^- 1 rv*^ \ 

= - [A + iBxix ■ h iB 2 j 2 . 671. — 72. — 7ll72l7i-72 

2 V r r r r J 

1 f 0C~^ 0C^~* 0C~^ 0C~^ ~\ 

+ -7il72l7i-72 A + iBffi.— + iB 2 7 2 T C^.—^.— . (3.16) 

2 V r r r r / 

The difficulty of the calculation stems from the fact that the matrices .x T and 7J .x T 
do not commute with 71.72. 

To proceed further, we introduce in the subspace of x T a longitudinal direction, parallel 
to x T , and a transverse plane, orthogonal to it. We define: 

T 

^ , x T2 = -1 , 



^ = 


T 


V-x T2 


T _ 
7a/i 


7a/i "I - 7a/Li 


7a£ = 


T ~T 
7a , 



r 



la-X X p + 7 a(U ; la- X ) 

7a 2 , = -1 (a = 1,2) . (3.17) 



The capital indices L and T [Eqs. (2~T3)] concern the longitudinal and transverse com- 



ponents with respect to the total momentum P, while the small indices I and t concern 
those of the three-dimensional relative distance x T . We list here some useful relations 
satisfied by these matrices: 

71-72 = 71L72L + ll-l 2 = 11L12L - lUl2i + ll-li , 

lUlll = -7lL72L7l5725^12 , ll-l 2 = -7lL72L7l5725^ , 
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7i-72 = -1u1m{w - wi 2 ) , 



'-fall, (W - W 12 ) 







(a =1,2). (3.18) 



([ , ] + is the anticommutator.) 

At the first stage of the calculation, one can eliminate 7 a 5 and 7* (a = 1, 2) in terms 
of 7 a £, 7 a £, w and Wi 2 . For the subspace of the matrices 7^, one introduces the following 
projectors: 

P++ = t(1 +i7w)(l + «72£) , P+- = ^(l + i7i«)(l-*72<) , 

P-+ = j(l-i7i«)(l + i72<) , P— = i(l-z 7 w)(l-i72^) , (3.19) 



which allows the decomposition of the 7^'s along the latter. Potential Vb [Eq. ( p.16 )] 
takes now the form: 



1il12lV = -(Ajuj2<e - C')7il72l + (A - Cjul2e)(l + liLl2Ll[-ll) 

+i(B llie + B 2l2£ )(l + jtL^Lllll) ■ (3.20) 



Notice that the second term in the right-hand side above commutes with the two others 
and therefore its exponential can be factorized and calculated independently. The first 
and third terms can be written in terms of the projectors (|3.19|) : 



-AiurYztYiL^L + CliLl2L + i(-Bi7 U + B 2 ^ 2( ,)(l + liLl2LlW 2 ) 
= V ++ {A + C) llLl2L + {B 1 + B 2 )(l+ llLl2L (w-w 12 )) 

{-A + C) llLl2L + (B 1 - B 2 )(l - llLl2L (w - w 12 )) 

+V-+ (-A + C)7il72L - (Bi - B 2 )(l - 7il72l(^ - w 12 )) 

(A + C) llLl2L - (Bi + B 2 )(l + llLl2L (w - w l2 )) 



(3.21) 



The exponential of this expression can be calculated by a series expansion. The pro- 
jectors V ++ , etc., commute with 7il72l(u; — Wi 2 ), but satisfy particular commutation 
rules with 71L72L. One factorizes the projectors V ++ , etc., on the left of the series. Each 
multiplicative factor of etc., can be resummed into exponential functions. At the 
end, one rexpresses 7^ and 7 2 £ in terms of Jil, 7 2 l, 715, 725 and the spin operators fl3.15|) 
and one introduces back the projectors V\ (i — 1, . . . , 4), (|3.14|) . One thus obtains: 
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2V B = e 2 llLl2L V 



1(1 + w 12 )e a ^+(f ++ + g+^V, + 1(1 - ^ 12 )e a -^(/__ - ^_)^ 
+ ^(1 - Wl2 )e a +^ (/++ + ^ ++ )P 2 + 1(1 + «;i 2 )e a -^ (/_ - <?__)P 2 
+ 1(1 - w 12 )e a ^+(f + - - g + -)V 3 + 1(1 + Wl2 )e a -T+(/_ + + ^ + )p 3 
+ 1(1 + ^ 12 )e a +^- (/ + _ - g+-)T, + 1(1 - ^ 12 )e a -^- (/_ + + g^)V, 



+ -^{lU-l2l) 



e a +^+h ++ Vt + e a +^-h ++ V 2 + e a +^h + _V 3 + e a +^/i + _P4 



(3.22) 



The definitions of the potential functions are the following: 



f rs = COsh y «2 + /?2 7 2 ; ^ 
, &7s 



sinh \/a% + /3^7 



2^2 



a 2 + /3^ T 



2 



sinh J a 2 + /3 2 7 2 



r, s = ± 



f 2 _ Q 2 _ ^2 



(3.23) 
(3.24) 



a± = 2(A±C) , P ± = 2{B 1 ±B 2 ) , 
7± = l±(w- w u ) . 

The exponential function e~^ B is obtained from Eq. ( |3.22| ) by the replacements a — > —a 
and (3 — ► — /3; is obtained by the replacements a — > a/2 and (3 — > /3/2, etc.. Also 
notice the commutation relations: 



la£l± = l^lai , 1alf±± = f± T lal , 
lai9±± = 9±Tlae , 7a^±± = ^±t7q 



a = 1,2 



(3.25) 



Equation (|3.22|) and the similar ones with different arguments allow us to project the 
wave equations ( |3.7f ) or ( |3.10| ) with the aid of the projectors Vi (i = 1,...,4), ( |3.14j ), 
appearing on the utmost right of the expressions, on the four-component wave functions 
(|3.12|) or (|3.13|) . One thus obtains coupled equations for the four components ipi or 
V*Bi = 1, - . - , 4) and eliminating three of them one reaches a final eigenvalue equation 
involving only one of the components. 
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4 Determination of the potentials 

Before proceeding to the reduction of the wave equations, we shall determine the expres- 
sions of the various potentials appearing in Eqs. ( |3.16| ) and ( |3.22| )-( |3.24| ) in terms of the 
elementary Coulomb potential. 

Potential V has the structure ( |2.13| ) and is obtained from Eq. (|2.14|) with the substi- 
tutions (|2.11|) . We shall, for the moment, not use the particular expression ( |2.15|) of the 
initial potential An, but rather present the calculations for the general case. We assume 
that A is expressible as a power series of the elementary Coulomb potential a/(2Pj,r): 



A, 



A (v) 

a: 



oo 
n=l 



2P L r 



(4.1) 
(4.2) 



The substitutions ( |2.11|) yield for the fine structure constant a the modification: 

a — > a' = a(l - ibi"/u)(l - ib^e) , (4.3) 
which can be expressed in terms of the projectors ( |3.19| ): 

1 - 6i)(l - h)V ++ + (1 - & a )(l + 6 2 )P+_ 



a: 



a 



+ (1 + 6x)(l - 6 2 )P_ + + (1 + 60(1 + 6 2 )P__ 



(4.4) 



Since v [Eq. (|4.2j )1 is proportional to a, its modification is similar to that given by Eq. 
(|4.4|) . The initial potential A thus undergoes the change: 



An 



A. 



oo 
n=l 



v++(i - hni - b 2 r + v+4i - h) n (i + b 2 ) 



+V-+{1 + 6i) n (l - b 2 ) n + V — (1 + h) n {l + b 2 ) n 
V ++ A (v(l - 60(1 - 6 2 )) + V + -A Q (v(l - 60(1 + 6 2 )) 
+V-+Ao(v(l + 60(1 - 62)) + V—A (v(l + 60(1 + 6 2 

v++V— + TV v_ + + p_+v;_ + V—V++ . 



(4.5) 



Reexpressing the projectors V ++ , etc., in terms of the matrices 7^ (a = 1,2), one obtains: 

A = hy^ + + y + _ + V ++ ] + hy^ + V_ + - - V ++ )i 7ie 



V- 



V+_ + Vl_]i 2 7 W 7 2<? . 



(4.6) 
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This expression should be identified with the combination (A + iB^u + iBi'yu ~ CjuJm) 
appearing in Eq. (|3.16|) on the left or the right of 7il72L7i-72 (the substitution ( |4.4| ) being 
done symmetrically with respect to this operator). We then obtain the identifications: 



and from Eqs. ( p. 24 



A 
B 1 
B 2 
C 



V, + V, + V. 



++J ) 



V- 



v. 



(4.7) 



a: H 
At 



V- + v ++ , 
V— - v ++ , 



V-+ + V+. , 



(4- 



For the particular case of Todorov's potential (|2.15|) , the expressions of the potentials 
V ++ , etc., are: 

1 



V- 

v + 



4 
1 



In l+4v(l- 6i)(l-&2 



In (l + 4v(l-6i)(l + 6 2 ) ) , 



-In (l+4w(l + & 1 )(l-6 2 )) 



V ++ = -In (l + 4v(l + 6i)(l + b 2 ] 



(4.9) 
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5 Reduction to a final eigenvalue equation 



In the absence of anomalous magnetic moments, the wave equations ( |3.7| ) or ( |3.10| ) can 
be reduced to a single Pauli-Schrodinger type equation for the component fa or ^bz 



22 1 . A similar reduction can also be undertaken here; however, due to the complexity 



of the new terms in the effective potential [Eq. ( |3.22|) 1, the reduction process is not as 
straightforward as before. The reason is that the components fa or ip Bi (i — 1, • • • ,4) 
[Eqs. ( |3.12| )-( |3.13| )1 do no longer have, in the general case, simple characterizations with 
the quantum numbers I (orbital angular momentum) and s (total spin). For instance, in 
the absence of anomalous magnetic moments, the component fa can be classified according 
to the quantum numbers I = j ± 1, s = 1 (j being the total angular momentum) and 
I = j. In the present case, this simple property is lost and such a classification will 
concern combinations of ip3 and fa. 

It turns out that the most convenient representation where the reduction process can 
be achieved is the "anti-Breit" representation defined with the wave function transforma- 
tions x — e^ B ip = e B ipB [cf. Eq. (|3.9|) 1. In this case the reduced wave function is a 
tractable combination of X3 an d Xi- We shall not, however, present here the reduced wave 
equation in the general cases of quantum numbers, the corresponding expression being 
still lengthy, but rather shall content ourselves with the simplest case of the j = quan- 
tum number, corresponding to the two sectors of pseudoscalar and scalar states. These 
are also the most sensitive sectors involved in zero-mass bound state problems in strong 
coupling regimes. 

Actually, for these sectors, the Breit representation (|3.9| ) is the simplest one and it is 
sufficient to project Eq. (|3.10p along the components ip Bi (i — 1, • • • ,4) [Eq. (|3.13|) 1. In 
these sectors, the operators wi 2 , (Wis + W2s)-x T and (Wis + W2s)-P T have the following 
quantum numbers: Wi 2 = 1, (Wis + W 2 s)-% T = 0, (Wis + W 2 s)-P T = 0. Equation ( p.!0| ), 
together with Eq. ( |3.22|) , then yields the following four coupled equations: 

P L e a+1+ (f++ + g ++ Wm - K - m 2 )^ B2 + (^-)(W XS - W 2S ).p T ^ B3 

- l -P L (^ r )(Wi S -W 2S ).x T e a +^-h + ^ BA = 0, 



Bl 



+ 1 -Pl(^)(Wis - W 2S ).x T e a -^+h^ B3 
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P L e a - 1+ {f-+ + g-+)^B3 - M^j B4 + (-— )(W 1S - W 2S ).p T ^ B1 

-\PL^ L ){W xs -W 2S ).x T e a -^-h_^ m = 0, 

+ \Pl{j^){W 1s - W 2S ).x T e a ^+h + ^ Bl = 0. 



(5.1) 



These equations allow the elimination of the three components ipBi, ipB2 and ipsi in 
terms of ipB3, which is a surviving component in the nonrelativistic limit. Defining 

\2 



D 2h + _ _+_-._ ( m i - m *f _-(a-7 + + a +7 _) (/+- ~ /?+-) 
e -2u = e -a +7 - ( /+ _ _ g+ )e -2h + ^ + ) 



(5.2) 
(5.3) 



and making the wave function transformation 



*pB3 = e u 4> 3 



(5.4) 



one obtains the following eigenvalue equation for 03, written, for simplicity, in the cm. 
frame: 



P 2 



3 a_7+ + a+7_ 



M 2 



4 (/+-.-</+_)(/_+ -</_+) 4 (/ + _-s + -) 2 P 



(m 2 — m 2 ) 2 



1 



4M 2 (/_ + _^ + )2 
M h. 



+ 



(m 2 — m 2 ,) 2 



l + /,2 + _) e -(a-7+ + a+7- 



-4ft 2 x 2 
+ (6^ 2 - 4S 2 
+4ft 2 x 2 



4P 2 (/+_-(/+_)(/-+-»- 
(mi — m 2 ) 2 



4ftr (/+_ - g + _) ATir (/_+ - #_+) 

, M (mi — m 2 ) 



u" + 



Akr (/+_ - #+_) 4ftr (/_ + - 

M / \ (mi — 777.2) / 



4^ I r(/ + _ -</+_) 



4fr 



i 3 = . 
(5.5) 



Here, the prime designates derivation with respect to r (= x ) 



2V 



f = 

J ^„2 ' 



dr 2 



(5.6) 



S is the total spin operator, S = s x + s 2 , S 2 = 2h 2 s, s = 0, 1; the other operators and 
functions are defined in Eqs. ( |3.23|) - (|3.24|) , ( |B3|) -( [5~3"1) . The eigenvalues of the matrices 
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7± in the sectors with j — are the following: 

J = , I = , s = : 7+ = 3 , 7- = -1 ; 

j = 0,£=l,s = l: 7 + = "I , 7- = 3 . (5.7) 

The sector with £ = 0, s = corresponds to the pseudoscalar states, while the sector with 
£ = 1, s = 1 corresponds to the scalar states. 
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6 Properties of the eigenvalue equation 



We study in this section two aspects of the eigenvalue equation ( |5.5| ) concerning, first, its 
short-distance singularities, and, second, its nonrelativistic limit. 

6.1 Short-distance singularities 

The question that arises here is whether the presence of the anomalous magnetic moments 
has any influence on the short- distance singularities of the effective potentials present in 
the eigenvalue equation. It was already clear from the expressions of the substitutions 
(|2.11|) that the form factors b a (r) (a = 1,2) should be bounded in modulus by 1 in order 
not to destabilize at finite distances the bound state system. A detailed analysis of the 
eigenvalue equation is however necessary to reach a more complete understanding of the 
role of the form factors near the origin. We shall limit our study to the case of Todorov's 
potential Q2.15|) . 



For a matter of comparison, we rewrite Eq. (|5.5|) in the case when the anomalous 
magnetic moments are absent. Here, we have B\ = B 2 = C = 0, A = A , a + = a_ = 2A Q , 



(3 + = (3_ = 0. Denoting h = |_ [Eq. (|5.2j)1 in this case, Eq. (|5.5|) becomes p2 



^l p 8Ao _ MlAA _ (ml -ml) 2 a Aq (m 2 - m 2 ) 2 
A 4 AM 2 ' AP 2 

-p 2 - Ah 2 * 2 !!' 2 + 6h 2 ti + Ah 2 * 2 h" 



-4S' 



(2A' Q + h')(l + 4x 2 A' ) - (A' + 2x 2 A' ') ] | 03 = 



(6.1) 



It is sufficient to study the short-distance singularity problem in the equal-mass case 
(mi = m 2 , h = 1) and in the ground state sector (£ = 0, s = 0). The dominant singularity 
comes from the term ^e 8A °, which, according to the expression ( |2.15| ) of Aq, yields the 



attractive potential a 2 /r 2 . This term is at the origin of the fall-to-the-center phenomenon 
with a critical value of a equal to | p4j. 

The above analysis can be repeated with Eq. ( |5.5| ). In the equal-mass case, one has 
bi = b 2 = b [Eqs. ( |2.11| )1 and the expressions of the various potentials [Eqs. ( p.23| )-( p^24|) , 
|)] become: 
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a+ = — In 
4 



1 + 4^(1 -6) 2 1 + 4^(1 + 6) 



n_ = - In ( l + 4u(l - 6 2% 



At 



l + 4v(l -6) 5 



l + 4w(l + 6) 2 



5+ _ = cosh ^/a 2 ^ + /^7 2 



0- 

Q!-l- 



ai + /?i 7 2 



sinh v/a 2 + /3+7 2 



—a. 



(6.2) 



Their behaviors near the origin are: 



2 l P L r' 



iln 14 

4 U + 6. 



< 6(0) < 1 



Ph 2 



3 a+ 



/+ _-, + _,e- + (l-^) + 

~ (PLr) 1 / 2 + ^ 7 !(ln(-^-))^(P L r)-^ , 
/-+-^- + ^(PLr) 1/2 . 



(6.3) 



The behavior of (/+_ — 9+-) near the origin crucially depends on that of 6(r). If 
6(0) 7^ 0, then /3 + (0) 7^ and hence (/ + _ — g+_) essentially behaves as r -1 / 2 . The first 
term in Eq. ( |5.5| ) has therefore a behavior of the type r _1 , contrary to the behavior 
of the type r~ 2 obtained in the absence of anomalous magnetic moments. Therefore, a 
non-vanishing of the form factors b(r) at the origin drastically changes the singularity 
of the effective potential at the origin. Also in this case, for s — 0, the function u [Eq. 
( |5.3| ) behaves as — a + and the combination — 4:h 2 x 2 u' 2 + 6h 2 u' + <ih 2 x 2 u" of Eq. ( |5.5| ) 
has a behavior close to h 2 /(Ar 2 ), which was absent in the initial case. This singularity is 
independent of the value of the coupling constant a and is located at the critical point. 
This would mean that the system, even for small values of a, would face strong attractive 
singularities, which are not observed experimentally. 

The above study suggests that the form factors 6(r) should vanish at the origin, in 
order not to drastically modify the situation found in the absence of anomalous magnetic 
moments. A smooth contribution of the anomalous magnetic moments would require 
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that the second term in the right-hand side of the equation of (/+_ — g+~), Eq. (p.'3\), be 
nondominant in front of the first. This implies that (3+, and hence b, vanish at least as 
rapidly as r 1 / 2 at the origin. (Also, in this case, the function u vanishes at the origin.) 

A parametrization of b a (r), corresponding to a vanishing at the origin as r, is obtained 
with the following choice of the functions c (r) of Eqs. (|2.12j) : 

Tiot 

c a (r) = d a + f a - (a = 1,2), (6.4) 

2irm a r 

with d a and f a constants, d a > 1, f a > 0. 



6.2 Nonrelativistic limit 

When the magnetic moment form factors b(r) are smooth functions, then for values of 
the coupling constant a of the order of 1/2 (the critical value), their effects can still be 
estimated perturbatively, their order of magnitude being fixed by a /it. An even cruder 
estimate is obtained by the nonrelativistic limit, which allows us to have easily an idea of 
the signs of the energy shifts. We shall assume that a is sufficiently small to also justify a 
perturbative-nonrelativistic treatment of the Coulomb potential v [Eq. ( [Op ] appearing 
in expressions concerning the anomalous magnetic moments. 

We treat the form factors b a (r) to first order. Among the effective potentials Q3.23 )- 



(|3.24 ) and (|4.8| )- (|4.9| ), only (3± and the h's are first order quantities in b a . In this approx- 



imation, the latter are given by the perturbation theory result: 



W'^Saib (a = 1 ' 2) - (6 - 5) 

One finds for the effective potentials: 

h+- - /3+7_ , h-+ ~ /3_7+ , 

13+ ~ -2v(h + b 2 ) , (3_ ~ -2v{h - b 2 ) . (6.6) 

Let the first-order perturbation due to the anomalous magnetic moments appearing 
in Eq. (|5.5| ) be represented by — 5V. We have: 



SV = (6h 2 - 4S 2 ) 
+4^ 2 x 2 



M _ K-m 2 ) 



M (mi - m 2 ) , 

-Al4— : 



Ahr Ahr 
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V 2 

-4S 2 



dr 2 

M 
Ahr 



M 
Ahr 
(mi 



m 2 



mi - m 2 ) 



Ahr 



Ahr 



■h. 



where V 2 is the laplacian operator. Using the expressions of the h's and fe's [Eqs. 
(|6^q) 1 we also have: 



(6.7) 

3)- 



M 
Ahr 



-h. 



(mi - m 2 ) 
Ahr 



h. 



Ahr 



M2v(b 1 + 6 2 )7_ + (mi - m 2 )2v{bi - b 2 )i^ 



1 



a: 



■ (^M 2 7_ + (mi 



m 2 ) 7+ 



(6.* 



167r^, mim 2 M r 3 

Then the corresponding perturbation in the nonrelativistic hamiltonian, designated by 
5Vjv_r, is related to SV with the relation [021 : 



5V 



2mim 2 
M 



5V, 



NR. 



(6.9) 



Introducing the total spin quantum number s (= 0, 1), we obtain: 

•• 2 1 



SV, 



a 



NR 



Ahm\m\ 
a 2 1 



^(x)[M 2 7_ + (m 1 -m 2 ) 2 7+] 



■[M 2 7_ + (mi -m 2 ) 2 7+] 



(6.10) 



Airh m\m\ r 3 

The first term contributes to the sector with s = 0, i = 0, for which 7 + = 3, 7_ = — 1 
[Eqs. Q5-7D 1, while the second one to the sector with s — 1, £ — 1, for which 7 + = — 1, 
7_ = 3. The energy shift then becomes: 



6E 



a mim 2 
2^ M 3 

a 5 mim 2 
~6^ M 3 



m 2 + m 2 — 4m!m 2 ) 



m 2 + m 2 + 4m!m 2 



SeoSso 
ni 



(6.11) 



with = £ + n' + 1, n' = 0, 1, . . .. 

While the sign of the energy shift is negative for the sector with s — 1, I — 1, it 
depends on the ratio mi/m 2 for the sector with s = 0, £ = 0. For the particular case of 
equal masses, m x = m 2 , the energy shift for the latter sector is negative and equal to the 
energy shift of the sector with s — 1, £ — 1. 

In the infinite mass limit, m 2 — > oo, the problem reduces to that of a spin-^ particle 
with anomalous magnetic moment placed in an external static Coulomb field. In this 



21 



case, the sector with s = 0, I = tends to the new sector with j = |, £ = and the sector 
with s = 1, I = 1 to the sector with j = ~, £ = 1. Equations ( |6.10| ) and ( 6 . become: 



5V/ 



ATfl 



2mf (X) ~ 2^f 5 V 



rrn a 5 a 5 <5j i/2 

o£/ = mi ^ mi- 



27r ra» 



67T n« 



(6.12) 



They agree, as they should, with the corresponding formulas obtained directly from the 
Dirac equation |25|. (Comparisons of theoretical predictions involving anomalous mag- 
netic moments with experimental data can be found in Ref. f26|.) 
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7 Summary and concluding remarks 



Using a matrix substitution rule, applied to the electric charge and deduced from the 
lowest order contribution of the vertex correction in QED, we introduced in local form 
the anomalous magnetic moments at each vertex of the higher order terms of the con- 
straint theory fermion-antifermion interaction potential. Since the latter already has a 
local form in three-dimensional rr-space, determined from summation of infra-red leading 
terms of multiphoton exchange diagrams, the new potential that arises also has a similar 
locality property and is calculated by a resummation of the corresponding series after the 
incorporation of the anomalous magnetic moments into the vertices. 

Focusing our attention to the sectors of pseudoscalar and scalar states (j = 0), the 
corresponding wave equations were reduced to a single eigenvalue equation. The re- 
quirement that the short-distance singularities of the effective potential should not be 
drastically enhanced by the presence of the anomalous magnetic moments imposed on 
the accompanying form factors the condition of a sufficiently rapid vanishing at the origin 
(faster than r 1 / 2 ). It is expected that when this condition is realized, then the incorpora- 
tion of the anomalous magnetic moments, even in the strong coupling regime, should not 
introduce qualitative or destabilizing changes in the properties of the fermion-antifermion 
bound state system. 
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